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The MILC collaboration's simulations with improved staggered quarks are being extended with 
runs at a lattice spacing of 0.06 fm with quark masses down to one tenth the strange quark mass. 
We give a brief introduction to these new simulations and the determination of the lattice spacing. 
Then we combine these new runs with older results to study the masses of the nucleon and the 
Q. in the continuum and chiral limits. 
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1. Introduction 

In this work we update our determination of the lattice spacing in the MILC collaboration's 
program of QCD simulations using A/j = 2 + 1 flavors of dynamical staggered quarks with the 
fl^^^ action[|l|], and update our calculations of the nucleon and Q. masses. Determining the lattice 
spacing is central to any simulation program. We emphasize that the same lattice spacing determi- 
nation, with the same uncertainties, affects everything computed from these lattices. The spectrum 
of the light quark hadrons is an important test of a lattice QCD simulation, and is a standard part 
of all large lattice QCD programs. Getting these well known masses right gives us confidence in 
our ability to predict unknown masses and matrix elements. Both the nucleon and Q. are stable 
against strong decays and well known experimentally, and hence are good tests of our program. 
Our previously published work on the light quark spectrum and lattice spacing determination is in 
Refs. [gi. 

We are now doing simulations with a smaller lattice spacing, a ?a 0.06 fm, in addition to 
lengthening some runs and adding runs at other quark masses at larger lattice spacing. The runs at 
a K, 0.06 fm are all still in progress, so all results reported here are preliminary. Table I shows the 
simulation parameters for the a 0.06 fm runs. 
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Table 1: Lattice parameters. Lattice spacings come from the smoothed ri/a values and ri = 0.318 fm. 
"RHMD-«" is a RHMD algorithm using n pseudofermion fields. Details of the algorithms will be given 
elsewhere. 



2. Finding the lattice spacing 

The results of a lattice computation are in units of the lattice spacing a, so a determination 
of a in physical units is essential. This requires taking some physical quantity (or average over 
quantities) as a standard. At this time we are using the T mass splittings determined by the HPQCD 
group[Q]. A possible alternative is fj^, and we are investigating whether this would be better[Q]. 
Because T mass splittings are not available on all of our ensembles and because the static quark 
potential can be accurately determined, we use ri defined by rjF(ri) = 1.0 0] as a quantity to 
interpolate the lattice spacing, as described in Ref. [^]. Thus, the T splittings give n = 0.318(7) 
fm , which is then used to determine the lattice spacing for all gauge couplings and quark masses. 

Static quark potentials were calculated from correlators of time direction lines, with Coulomb 
gauge fixing for the spatial connections!^. As the lattice spacing decreases, determining the static 
potential becomes more difficult because of the increasing constant term, corresponding to the sin- 
gular self-energy of the static quark. Since the Wilson loop expectation value falls off exponentially 
as {W{R,T)) f» exp (— y(/?)r), a larger constant in V{R) translates into an exponentially smaller 
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expectation value of lV(/?,r), and larger fractional statistical error. This can be alleviated, at the 
cost of some spatial resolution, by smearing the time direction links in the lattice. 

Figure la shows the static quark potential, with and without one iteration of APE smearing 
of the time direction links, while Fig. lb shows the static quark potential at mq = QAms for four 
different lattice spacings. 
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Figure la: The static quark potential, with 
and without APE smearing of the time di- 
rection links. This is for the aniq/anis = 
0.0072/0.018 run. The constants in the po- 
tential differ (that's the point), so we added 
0.27731 to the smeared potential to show the 
agreement. (This shift forces the potentials to 
agree at ri.) 
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Figure lb: The static quark potential at 
niq = OArris for four different lattice spacings, 
roughly a = 0.06, 0.09, 0.125 and 0.15 fm. 
The a = 0.06 fm potential uses smearing of 
the time links. The colored rulers near V = 
are in units of the corresponding lattice spac- 
ings, while the "physical" ruler near V = —2 
is in units of 0. 1 fm. 



In Fig. ^ note the visible lattice artifact at the on-axis (2,0,0) point marked with |!. The 
spatial range used for the r\ determination at a f» 0.06 fm was 4.01 to 7.01. Also note that the 
a = 0.06 fm potential (red) at r = (marked with ^!) is much smaller in magnitude than if it 
continued the trend, again showing the effect of the smearing. 

To find ri the potential was fit to y(r) = C + A/r + Sr for a range around r\. With a 0.06 
fm, lattice artifacts are unimportant in these distance ranges. Errors are from a jackknife analysis, 
using a block size of 20 lattices, or 100 molecular dynamics time units. Figure 2a shows estimates 
of r\ I a from various spatial and temporal distance ranges without smearing, while Fig. 2b (right 
half) shows fits with one iteration of APE smearing, including projection onto SU(3). The fits 
shown in Figs. 2a and 2b (right) are from fitting the "Wilson loops" to a single exponential: V(r) = 
log(lV(r, r)/W (r, r + 1)). We also do fits including an excited state of the potential. This gives less 
dependence on fit range (smaller systematic error) but larger statistical errors. The r\ja estimates 
on the left side of Fig. 2b are from such two-state fits. Fitting this, together with our previous 
results at a PS 0.018, 0.015, 0.0125 and 0.09 fm, we get a "smoothed r\ \ which is then used to set 
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Figure 2a: Estimates of ryja from "loops" 
with unsmeared time links. This is for the 
amq/anis = 0.0072/0.018 run. 
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Figure 2b: Estimates of r\/a from the same 
run, with one iteration of APE smearing in- 
cluding projection onto SU(3) for the time di- 
rection links (right side). Fits including an ex- 
cited state (left side). The arrow indicates the 
point taken as our estimate. 



the length scale for all ensembles. 

3. Nucleoli and Q. masses 

Both the nucleon and the Qr are stable against strong decays, and so their masses should be 
accessible in a straightforward way to lattice calculations. In principle the nucleon mass could 
be used to determine the lattice spacing. In practice its correlators are noisy, and the non-trivial 
chiral extrapolation makes an accurate determination difficult. The Q. has a very mild chiral 
extrapolation — no terms with log(m;i:) or m\ at one loop order. Since it contains three valence 
strange quarks, its mass is very sensitive to the strange quark mass. Therefore, the Q. mass checks 
our determination of the lattice strange quark mass from pseudoscalar physics. 

Figure 3a shows the nucleon correlators from the a « 0.06 fm runs at 0.4/?!^ and Q.lnis. Note 
the expected difference in the slope, the alternating (opposite parity) signal at short distances, 
and the increasing error bars at large distances. Figure 3b shows fits to one of these correlators 
{niq K, Q.lrris) as well as correlators for roughly the same light quark mass aia^ 0.125 and 0.09 fm. 
In this plot the size of the symbol is proportional to the confidence level of the fit, with the symbol 
size used in the captions corresponding to 50%. From graphs like this we look for a plateau and 
reasonable confidence level {x^). The arrows in the graph indicate the fits that were chosen. Most 
of the difference in the masses comes from the fact that the light quark masses in these three runs 
were not exactly equivalent, and some comes from effects of lattice discreteness. 

In Figure 4a we show nucleon masses for a ^ 0.0125 fm, 0.09 fm and 0.06 fm, where the three 
colored arrows in Fig. 3b indicdate the fits that were chosen to include in this figure. The black 
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error bar at the left is the experimental value. The error on this comes from the uncertainty in r\ 
(We use 0.318(7) fm). Another point at a « 0.06 fm at m/ ?a O.lmv ((m;r''i)^ ~ 0.12) is in progress. 
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Figure 3a: Nucleon correlators at a = 0.06 fm 
for light quark masses 0.4 /m^ and 0.2 m^ . 
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Figure 3b: Nucleon mass fits for m; = 0.2ms. 
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Figure 4a: Nucleon masses (red blue and 
green symbols), and fits to the nucleon mass 
(cyan and magenta). 
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Figure 4b: The same, expanded to show the 
small quark mass region. The points with 
mi = O.lnis are from the fits indicated with 
arrows in Fig. 3b. 



We show two forms for extrapolating the nucleon mass to the physical quark mass and contin- 
uum limit. The first (cyan) includes only the lowest order nonanalytic correction, order m^. 



rriNn =Mori+Ai{m„ri)^+Bia^a + r, ^ m\ 



(3.1) 
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In these fits we just set fji, gA and Ma — Mp^ to their physical values. The cyan and magenta 
symbols in Figs. 4a and 4b are the data points with the fitted lattice spacing corrections, Bio^a, 
subtracted, i.e. at a = 0, and the cyan and magenta lines are the fitting function at a = 0. 

The second form (magenta) [^] includes the A without the assumption that A = M/^ — M^ » 
fji. When Taylor expanded in powers of mji this includes terms up to m^, m^log(m;t). 

rriNn = Mon +Ax{mnnf +A2{mnnY +Bia^a +B2{mjcnfa^a + n :^^j^H{mjc,Mt, -Mn,A) 

(3.2) 

where, if m^;: < A = Ma — 

H{m^,A,A)=Ahogi2A/m^) + Aml(^^log{m^/A)-l^ (3.3) 
- (a2 - m2)'/'log (^A/rn, + ^A^ml - l) - Am^ 01og(2A/A) - ^ 
if m;t > A 

//(m;„A,A) =A3log(2A/m;,)+Am2 Qlog(m;,/A)-l^ (3.4) 
- (m2 -A2)^/^arccos(A/m;,) -Am2 01og(2A/A) - 

Systematic errors have not been properly analyzed yet, and the data set is incomplete. With 
these reservations, the order chiral fit in Fig. 4 at the physical quark mass is M^ri = 1.52 it 
.04 lb syst., or MN^MeV) = 942 it 25 ± syst., where systematic errors on M^ri must be at least 
0.045, the difference between the two chiral fitting forms. We are also experimenting with chiral 
extrapolations using fitting forms from partially quenched staggered chiral perturbation theory [^]. 

In fitting the D. mass the chiral extrapolation is simpler, but we must deal with its dependence 
on the strange quark mass. Our lattice simulations were run with an estimate of the strange quark 
mass and only after running the simulation cam we determine what the correct strange quark mass 
is. Thus, we compute correlators at two valence strange quark masses, and then interpolate to the 
valence strange quark mass determined from pseudoscalar mesons (basically M/f). Figure 5a (like 
figure 3b for the nucleon) shows fits to the D. correlators for runs at three different lattice spacings 
with nil ~ 0.2ms, again with arrows indicating the fit we chose. For the chiral and continuum 
extrapolation we use a linear extrapolation in ami and a^a, which fits well with /dof = 6.0/1 . 

MQ,ry=MQry+A{mj,rif+Ba^a (3.5) 

Figure 5b shows the Q. masses interpolated to the strange quark mass determined from pseu- 
doscalar mesons as a function of light quark mass. It also shows the fit form above for a 0.125, 
0.09 and 0.06 fm, and in the continuum limit, i.e. with B set to zero. The black symbol is the 
experimental value for mQr\ , where the error is from the uncertainty in ri . 

Errors on the extrapolated lattice result are statistical (0.7%), determination of from pseu- 
doscalars (0.6%), and an estimate of neglected NNLO terms in chiral extrapolation (qos)- Other 
small errors come from using the wrong strange sea quark mass, and the choice of procedure for 
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fixing strange quark mass. Using data available at conference time, the continuum and chiral ex- 
trapolated value is M^j * ri = Converting to physical units using r\ = 0.318(7) and 
adding errors, this is Ma{MeV) = 1660 ^^q, where the experimental value is 1672 MeV. 

This work is supported by the US Department of Energy and National Science Foundation. 
Computations were done at the NSF Teragrid, NERSC, USQCD centers and computer centers at 
the Universities of Arizona, Indiana, Utah and California at Santa Barbara. 
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Figure 5a: Q. mass fits for mi = 0.2/Ms. In 
these plots the strange quark mass has not 
been tuned to the correct ms', it is one of the 
masses at which correlators were computed. 
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Figure 5b: Q. masses and fitting (see text). 
The error on the experimental point is from 
the error in ri . 
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